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Abstract 

Despite the recent evidence that anti-de Sitter spacetime is nonlinearly unstable, we argue 
that many asymptotically anti-de Sitter solutions are nonlinearly stable. This includes 
geons, boson stars, and black holes. As part of our argument, we calculate the frequencies 
of long-lived gravitational quasinormal modes of AdS black holes in various dimensions. 
We also discuss a new class of asymptotically anti-de Sitter solutions describing noncoa- 
lescing black hole binaries. 
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1 Introduction 

It has long been known that anti-de Sitter (AdS) spacetime is stable to linearized perturbations, 
and positive energy theorems ensure that AdS cannot decay. However, there is growing evidence 
that AdS is nonlinearly unstable, in the sense that any finite perturbation will eventually 
become large and, at least in many cases, form a small black hole [U [21 E]- There were at least 
two earlier general arguments for this instability, which take the following form. Physically, 
AdS is like a confining box. If Einstein's equations are sufficiently ergodic, then one might 
expect any generic finite perturbation to eventually explore all configurations consistent with 
the conserved quantities, including small black holes [I]. A more mathematical argument is 
that the linearized perturbations do not decay. So nonlinear corrections (which usually involve 
integrals of powers of the linearized modes) are expected to become large [3]. Here, as is 
standard in the relativity community, our use of the term AdS refers to the universal covering 
space of the surface Yli=i{X 1 ) 2 — Tf — T| = — L 2 in a Minkowski space of signature (d — 1, 2). 
This is what string theorists call "global" AdS, as opposed to just the Poincare patch of AdS 
which is often used in studies of the AdS/CFT correspondence. 

The recent evidence for the nonlinear instability of AdS came first from numerical evolution 
of a spherical scalar field coupled to gravity [U [2|. It was found that no matter how small 
one makes the amplitude of the scalar field, the energy in the scalar field cascades to smaller 
and smaller scales and eventually forms a black hole. The time to form the black hole tsH 
was found to scale with the initial amplitude e via tsH ~ V e2 - The instability (and this time 
scale) can also be seen in terms of resonance^] in the perturbative construction of the solution 
PP. For a generic initial perturbation, resonances appear at third order which require that one 
add higher frequency modes with linearly growing amplitude. In other words, the perturbative 
solutions cf) have the schematic form 

4> = <f>o + A4>l + tB4>l + ... (l.i) 

Resonances arise when there are nontrivial linear relations among the frequencies of the first order equations. 
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where 0o is a solution to the linearized problem and A, B are appropriate operator^} For 
t ~ l/0o , the term linear in t becomes comparable to the lowest order term. After this time 
one expects nonlinearities to dominate. Since total energy is conserved, the energy is transferred 
to the higher frequency mode. Here it is important to note that, in describing the above result, 
we have already included the possibility to resum perturbation theory by shifting the mode 



frequencies at each order. The expansion takes the form (1.1) even after all possible such shifts 
have been implemented. One expects the transfer of energy to higher frequencies to continue, 
eventually forming a black hole. 

The purely gravitational problem has a similar perturbative structure [3J. For a generic 
initial perturbation, resonances again arise at third order which transfer energy to higher fre- 
quency modes. Although the numerical evolution has not yet been performed for this case, it 
is natural to expect eventual black hole formation in this context as well. 

Although generic perturbations of AdS lead to a nonlinear instability, it was shown in [3] 
that if one starts with certain individual graviton modes, nonlinear corrections can be added 
without triggering an instability. The resulting nonlinear analog of the linearized mode is called 
either a georj^] or an oscillon depending on whether or not it has a Killing field that is timelike 
in some region. Given the linearized mode, there is a one parameter family of either geons or 
oscillons labelled by their energy. They can be viewed as gravitational analogs of boson stars 

The two general arguments for the nonlinear instability of AdS mentioned above would seem 
to apply to geons and boson stars as well, suggesting they too will be nonlinearly unstable. We 
will argue that this is incorrect: Most asymptotically AdS solutions which are stable to linear 
perturbations appear to be stable nonlinearly as well. Strictly speaking, the lack of a useful 
Killing field means that our methods do not directly address the stability of oscillons. However, 
we expect their behavior to be similar to that of geons and other solutions with some notion 
of time translation symmetry. 

The key difference between pure AdS and these asymptotically AdS solutions is that in 
AdS, the normal mode frequencies are all integer multiples of the AdS frequency. This results 
in a large number of resonances. Other solutions with less symmetry do not have this prop- 
erty. However, they still have normal modes with very large angular momentum i which are 
supported at large radius due to the centrifugal barrier. In this asymptotic region, the met- 
ric approaches AdS and the normal mode frequencies approach those of AdS. In other words 
there are approximate resonances at large I. Indeed, the purely gravitational instability of AdS 
seems to rely on these large I modes. Since a generic perturbation always contains some large 
i modes, one might wonder if the AdS instability carries over to these solutions as well. 

By studying the perturbative expansion in more detail, we will argue that these approxi- 
mate resonances are not strong enough to trigger an instability In particular, geons appear to 
be nonlinearly stable. Since geons are continuously connected to pure AdS by decreasing their 
amplitude to zero, one obtains an interesting picture of stability in the space of asymptotically 
AdS solutions. The geons radiate out from AdS, and are stable in finite neighborhoods whose 
size decreases to zero as one approaches AdS (see Fig. [T|. In other words, there are "islands 
of stability" in a neighborhood of AdS. A similar picture of stability has recently been rigor- 
ously established for a simpler nonlinear PDE on a compact space: The nonlinear Schrodinger 
equation on a torus [HJ [91 [TO] . 



2 For example, the term A<\>\ represents a a trilinear operator with 'diagonal' arguments </>o, </>0j </>o- 
3 This use of the term geon should should not be confused with the use of the term in, e.g., [6]. 
4 More precisely, we will argue that instabilities arise only in high dimensions for perturbations of low differ- 
entiability 
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Figure 1: The structure of solutions near empty AdS. The shaded regions denote islands of 
stability near 1-parameter families of geons or oscillons, both indicated by black lines with 
arrows pointing toward increasing amplitude. Because perturbation theory about empty AdS 
leads to geons only for a measure zero set of seed solutions, each such region has been drawn 
so that empty AdS lies at a cusp. 

Given the stability of general solutions without horizons, AdS black holes that are linearly 
stable appear to be nonlinearly stable as well since the modes eventually decay (although the 
large i modes live a very long time). However, the AdS boundary conditions do give rise 
to a novel class of black hole solutions: noncoalescing binaries. We will argue that there 
are asymptotically AdS solutions describing two black holes in a binary orbit which does not 
decay. The radiation produced by the black holes bounces off infinity and stops the orbit 
from shrinking, i.e. the gravitational radiation is in a standing wave configuration. Although 
these noncoalescing orbits can be stable against radial perturbations, at least one black hole is 
unstable to superradiance. So eventually, the angular velocity is reduced and the black holes 
spiral together. Similar orbits around the Kerr black hole were proposed in [11], and recently 
realized in [12], even though they are based on a different mechanism. 

We begin in the next section with a brief discussion of the nonlinear Schrodinger equation 
on a torus. We will argue that the stability of asymptotically AdS solutions is similar to the 
structure found in this simpler system. To discuss the role of approximate resonances, we need 
to know the rate at which the normal modes with large angular momentum (or quasinormal 
modes for AdS black holes) approach those of AdS. This is computed in section 3 using a WKB 
approximation. We consider arbitrary dimension and both rotating and nonrotating black holes. 
Section 4 contains the main argument that perturbation theory converges despite the approx- 
imate resonances. Section 5 describes the construction of noncoalescing binaries. The final 
section discusses some generalizations of these results and their implications for gauge/gravity 
duality. 

2 Nonlinear Schrodinger equation on a torus 

Consider the cubic nonlinear Schrodinger equation on a torus of dimension D > 1 

id t u + V 2 u = X\u\ 2 u. (2.1) 
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This equation has a conserved energy 



ll> . I 1 IV7„.|2 i ^ i.i 1 



E = J d u x [^u\ z + -M 4 J (2.2) 

and a conserved norm N = J d D x\u\ 2 . For A > 0, solutions do not blow up in finite time. Given 
smooth initial data u(x,t = 0), there is a smooth solution for all t > 0. 

It has recently been shown that the solution u = is unstable in the sense that there are 
nearby solutions in which the energy cascades to higher and higher Fourier mode^] [HI E] • This 
can be made precise by looking at higher Sobelev norms defined as follows. Let u n (t) be the 
n th Fourier mode of u(x,t) on T D and define 

\Ht)\\%'= £(l + ny) s K(t)| 2 . (2.3) 

where rij are integers labeling the modes. It was shown that for any small number S and large 
number K, there is a smooth solution for which 11^(0)11^ < 5 and \\u(T)\\ 2 HS > K at some 
time T. This holds for any s > 1. (The s = 1 norm is constant due to energy conservation). 

This is directly analogous to the instability of AdS with the one exception that there are 
no analogs of black holes in this theory to cut off the cascade. The reason for the instability is 
also the same. Solutions of the linear equation (A = 0) are of the form 

u(x,t) = Ae i{n ^-^ (2.4) 

with oj = njV? . Thus there are lots of resonances. 



A simple class of nonlinear solutions to (2.1) are plane waves. These take exactly the 



same form as the linear solutions (2.4) but now the frequency depends on the amplitude, 
u> = rijinP + XA 2 . They are analogs of the geons in the gravitational problem. Faou, Gauckler 
and Lubich have recently shown that for almost all A, these nonlinear solutions are stable 
and there is no energy cascade JTU]. More precisely, for s > D/2, the H s Sobolev norm 
remains small (at least for times that extend to arbitrary negative powers of the size of the 
perturbation). The proof relies on the fact that for almost all choices of A, most frequencies 
are now nonresonant. There remain a finite number of resonances, but they can only transfer 
energy among themselves. Since the energy cannot leak to high frequency modes, the higher 
Sobolev norms do not grow. 

The picture that emerges is that if a perturbation of u = has one mode with amplitude 
much larger than all the others, then the solution will stay close to the plane wave associated 
with that mode. Otherwise the modes compete and nonlinearities drive the energy into higher 
and higher frequency modes. We will argue that stability issues in AdS are similar. 



3 Asymptotic behavior of quasinormal modes 

To study the nonlinear stability of general stationary backgrounds such as geons^j boson stars, 
and black holes, we need to know how quickly the frequency of the large angular momentum 
perturbations approach those of AdS. Since these modes are located at large radius, we may 
expect the leading deviation from the AdS frequencies to be determined by parameters like the 
total mass M and the total angular momentum J of the background. While the detailed form 



5 This doesn't happen in one spatial dimension since (2.1) is integrable in that case. 
6 We include geons here even though their Killing field is spacelike at infinity. 
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of these frequencies may also depend on higher moments of the energy-momentum distribution 
(which contribute to the metric at the same order in AdS), we may expect the discrepancies 
Au n = u n — u^ dS from the AdS values to decay at large angular momentum £ with the same 
power law for all values of these parameters^} This should be a consequence of the well-known 
Fefferman- Graham expansion for asymptotically AdS spacetimes^} 

We test this hypothesis (and determine the associated power law) in two steps below, first 



using the global AdS^ Schwar zsch ild black hole in section 3.1 , and then using a simple rotating 

d-3 



black hole solution in section 



3.2 



In both cases we find Re Au n ~ Kd£~~2~ at large £ where 
only Kd depends on M and J. This analysis also displays the expected long-lived quasinormal 
modes first observed in the case of a massive scalar field in [131 [Hj. For completeness, we also 
treat such scalar field perturbations below. 

3.1 Long-lived quasinormal modes: Schwarzschild- AdS black hole 

We first briefly review the Kodama-Ishibashi harmonic decomposition of gravitational pertur- 
bations in the Schwarzschild- AdS black hole and discuss carefully the boundary conditions that 
preserve the global AdS asymptotics of the spacetime. Then we find the quasinormal mode 
spectrum of modes with large angular momentum using a WKB approximation. 

3.1.1 Perturbations of Schwarzschild- AdS and their boundary conditions 

Consider the Schwarzschild-AdS^ black hole, with horizon topology S d ~ 2 , that asymptotes to 
global AdSd, 

ds 2 = _ fdt 2 + _ +r 2 dn 2_^ with /=1 + ___™_ (31) 

where L is the AdS radius, r m is the mass parameter, dVL 2 d _ 2 is the line element of a unit radius 
(d — 2)-sphere and the horizon is located at r = r + (the largest real root of /). 

We want to study gravitational perturbations in this background that are constrained to 
obey the linearized Einstein equations subject to appropriate boundary conditions at the hori- 
zon and at the asymptotic boundary. Since our d-dimensional background is locally the product 
of a 2-dimensional orbit spacetime (parametrized by the {t, r} coordinates) and the base space 
S d ~ 2 , we can decompose perturbations h a b according to how they transform under coordinate 
transformations on the sphere S d ~ 2 . This harmonic decomposition was worked out in detail 



by Kodama and Ishibashi [To] . The most general perturbation of (3.1) can be decomposed 
into a superposition of three classes of modes: tensor, vector and scalar. This decomposition 
is particularly useful since it allows one to reduce the linearized Einstein equation associated 
to perturbations of (3.1) to a set of three decoupled gauge invariant Kodama-Ishibashi (KI) 
master equations for the KI master fields $^ which can be written in a compact form as, 9 

(□ 2 -^$(f(t,r) = 0, where j = {T,V,S}, (3.2) 



The imaginary part vanishes for any linearly stable horizon-free background. 

8 Backgrounds which differ from AdS at large r by terms larger than those associated with M, J, etc. should 
have Aw n ~ £ ~ a with smaller a. Since the approximate resonances are even weaker, they should be even more 
stable. 

9 Note that tensor perturbations exist only for d > 5 which we henceforward assume to be the case when 
discussing this sector. 
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for tensor, vector, and scalar perturbations, respectively. Here CI2 is the d'Alembertian operator 
in the 2-dimensional orbit spacetime, and the potentials Uj are further discussed below. They 
depend on the properties of the background, namely on the mass parameter r m and cosmological 
length L, and on the eigenvalues of the associated (regular) tensor, vector, and scalar harmonics. 
The latter are, respectively, given by 

A . = £.(£. + d _ 3 )_ c . ) where {(£ j ,c j )} = {(e T ,2),(£ v ,l),(ts,0)} (3.3) 

with integer £j obeying £t > 2, £y > 1, and £3 > 0. 

We also study perturbations of an axisymmetric massive scalar field \1/ in the Schwarzschild- 
AdSrf black hole background. This field obeys the Klein-Gordon equation, 

- /x 2 ) *(i, r, x) = , *(t,r,x)=r-" *%\t, r) S^ (x) , (3.4) 
where SVq(x) are the scalar spherical harmonics on the (d — 2)— sphere with eigenvalues Xs 



defined in (3.3). Note that the spectrum of £^ coincides with that of £3 and, because the 
background space-time is spherically symmetric, we can restrict our attention to axisymmetric 
modes without loss of generality. 



It is useful to introduce the tortoise coordinate 

"°° dr 



,. . < z < 00 , (3.5) 

that increases monotonically as we move from the asymptotic boundary, where z ~ L 2 /r — > 0, 
to the horizon where z ~ — ^7^7^ +00. Since our background is time-translation invariant, 

the fields can be Fourier decomposed in time as ^ (t, r) = e -iWj '*<I> t 



3 ' *-3 



In these conditions, the angular and time dependence of our gravitational and scalar field 



perturbations factors out, and studying equation (3.2) and (3.4) boils down to solving the radial 
Schrodinger equation, 

<9, 2 <U + K - U s ) = , with j = {T, V, S, tl>}, (3.6) 

for tensor (T), vector (V), and scalar (S) gravitational perturbations and scalar field perturba- 
tion (ip). The exact expression for the potentials {Ut,Uv,Us,U^} can be found in equations 
(5.6), (5.15), (3.2)-(3.8) of KI PS], and (3.6) of [13], respectively. 

Solutions of this equation have to obey appropriate boundary conditions. At the horizon, 
only ingoing modes are allowed. To enforce this condition, we change to ingoing Eddington- 
Finkelstein coordinates (appropriate to extending the analysis through the horizon) and de- 
mand regularity in this coordinate system. The metric of the Schwarzschild-AdS^ black hole 
asymptotically approaches that of global AdS^; we would like to restrict to perturbations that 
preserve this behavior. The KI master variables have the asymptotic expansion, 

$ 2UL^Q ~ $0 r ^ +aj + $ i r~^r~ a i , with a T = 2, a v = 1, a s = . (3.7) 

The linear differential map = ($^) that reconstructs the metric perturbations (in 
a given gauge) can be found in [15J. The requirement that these metric perturbations are 
asymptotically global AdS^ in the sense of [161 E]> imposes the conditions [18j [3]: 

$1 = - m $ , if j = S in d = 4 ; $ = 0, otherwise . (3.8) 
Ac — 2 
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These are also natural boundary conditions in the context of AdS/CFT: they allow a non-zero 
expectation value for the CFT stress-energy tensor while keeping fixed the boundary metric. 
Other boundary conditions that might be called asymptotically globally AdS were studied in 
[19J, but turn out to lead to ghosts (modes with negative kinetic energy) and thus make the 
energy unbounded below [2"U] . 

In the case of the massive scalar field, the asymptotic expansion reads [2TJ [22] 



2 ~(tp) i *o , $1 \ d-l , I(d-l) 2 



r 2 ^V,^L +0O ~^ + ^> where A ± = — ^ = y + /i 2 L 2 . 

(3.9) 

Stability of the AdS background (i.e. scalar energy bounded below) requires the mass of the 
scalar field to obey the Breitenlohner-Freedman (BF) bound /i 2 > A i2 | BF , where fJ?\ BF = 
—(d — 1) 2 /(4L 2 ). Furthermore, to have a well-defined initial value problem we must impose 
appropriate boundary conditions at r — > oo. For scalars with mass in the range ^ 2 \ BF < fi 2 < 
fJ?\ BF + l/£ 2 , there is a choice of boundary conditions that respect the full AdS isometries: 
one can impose either that decays as r~ A +, or that it decays as r~ A - (in this case with 
no term at order r _A+ ). According to the AdS/CFT correspondence, this choice dictates 
whether the operator dual to $^ has dimension A + or A_. We will consider both cases. For 
fx 2 > fx 2 \ BF + 1/ L 2 , only fast fall-off boundary condition $ = (9(r~ A+ ) leads to a theory with 
energies bounded below; the alternate boundary condition leads to ghosts with negative kinetic 
energies [20J. We thus choose the boundary condition $ = in this case. 



3.1.2 WKB analysis of quasinormal modes 

We are interested in the large angular momentum WKB limit, Ej — > oo, of the radial Schroedinger 
equation (3.6). Define 



LUj = pj Wj , Pj = Ej + 



d-3 



(3.10) 



so that (3.6) takes the form 
1 



V 



S"(z)+S'(z) 2 -[V-w 2 ] 



where V 



f 



(3.11) 



Note that V is independent of the perturbation sector we look at, whereas Xj{ z ) depends on the 
particular type of perturbation. We can solve this equation in a WKB approximation, which 
we review below, with expansion parameter pj 1 1. Introducing the wavefunction expansion 

Swj,pj — Yl ( k=oPj k Sj,k we find the leading order terms, 



Sj,o(z) = ±i 



w 2 — V(x)dx, Sj : i(z) 



1 



In [w 2 - V{i 



(3.12) 



± 



V"{x) 



[w 2 - V(x)] 



3/2 + 4 



V'(x) 



4 Xi(^) 



[w 2 - V(x)] 



5/2 



w 2 — V(x) 



dx . 



These leading terms are enough to illustrate a relevant property of this WKB system: the 
source term Xji z ) only contributes at second or higher order of the expansion, i.e. to Sj t k(z) 
with k > 2. For our purposes it will be enough to include only the next-to-leading order 
correction, i.e. k < 1. Thus we only need information concerning the potential V(z). A typical 
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Figure 2: WKB potential V(z) for r + = L and WKB allowed (J, III) and forbidden (77) 
regions for the corresponding particle orbits. 

behaviour for this potential is shown in Figure [2j At first sight one might worry that all 
information about the specific sector of perturbations is lost at this WKB leading order. This 
is however not the case: this information will be encoded in the asymptotic boundary condition; 



see (3.13) below. 

We are interested in long-lived quasinormal modes. These are bound states with frequencies 
L~ 2 < Wj < V(z M ) related in the WKB limit to particle orbits trapped in the exterior region 
limited by the asymptotic AdS boundary (z = 0) and the potential barrier at z = za (see Figure 
[2]). However, as solutions to a wave equation they can slowly leak through the potential barrier 
to the black hole horizon at z — > 00. The role of the large angular momentum I is therefore 
to provide the potential barrier that, as we will find, traps the waves for sufficiently long time 
to make these quasinormal modes long-lived before being absorbed by the horizon. In a WKB 
analysis we have to consider three regions: the allowed region I (0 < z < Za), the forbidden 
region II [za < z < z B ), and a second allowed region III (zb < z < 00), where Za, zb are the 
turning points defined in Figure [2] Here the terms allowed, forbidden, and turning point refer 
to the particle orbits associated with the given modes in the WKB limit. 

We want to find the quasinormal modes frequencies in the limit of large Pj using a standard 
WKB approximation. We start at the asymptotic boundary and move downwards toward the 



horizon. The potential Uj in (3.6) typically contains terms proportional to p 2 but also terms 
proportional to r 2 (that are inherited from the r 2 /L 2 factor in the metric function /). At the 
conformal boundary, z — » 0, these terms compete with each other and both contribute to the 
asymptotic solution in region I. Our task is to find the associated asymptotic solution of the 



Schroedinger equation that obeys the asymptotic boundary condition (3.8). More concretely, 
the asymptotic potential and solution obeying the desired boundary condition are (as observed 
earlier, note that these do depend on the particular sector of perturbations through Uj) 



tfiLo ~T 2 + W - X) 477 ' l-o - ^ Ju, [ Pj z^Jw]--), (3.13) 




where J v (x) is a Bessel function of the first kind and we defined 



d-1 



d-3 



d-5 



(d-1) 



(3.14) 



The solution for pj — > oo with finite z 7^ can be obtained by expanding (3.13) for large pjz 



while consistently promoting z 
the desired WKB solution in region / 



V\z^q fo dx\/Wj — V(x). This gives 



cos 



Pj 



w 2 j -Vdx-j{l + 2v j 



Defining 



and 



Vj = Pj 



~A 



71 



1 + 2v 



3) 1 



Qj(z)=w\ )-V{z) 



(3.15) 



(3.16) 



(3.17) 



we can rewrite (3.15) in the standard WKB form for region / in which the corresponding 
particle orbits are allowed, 



1/2^1/4 exp 

p/ Qj 



r Qj dx 



1/2^1/4 eX P 

Pj Qj 



-1 Pi 



/ t^j dx 



ZA 



(3.18) 



Note that (3.17) encodes the asymptotic boundary condition information. Since region 77 
(za < z < zb) is a forbidden region for the corresponding particle orbits, the WKB solution is 



a 



U) 



V 



exp 



-Pj 



\Qj\ dx 



+ 



a 



p) f2 \Qj\ 1/A 



exp 



Pj 



2/1 



\Qj \ dx 



(3.19) 



where the respective amplitudes are fixed by the standard WKB connection formulas around 
the turning point z = za, 



C[ j) = I (c { 2 j) e -^/ 4 + C^ e 1 ^ 4 



dP=i(cf? ) e- iir / 4 -CPe i */ 4 ) . 



(3.20) 



Finally, region 77/ (zb < z < 00) is again an allowed region where the horizon is encountered 
at z — > 00. Its WKB wavefunction is 



$0') 



in 



Pi 



+ 



pfQf 



exp 



-i Pj 



(3.21) 



The WKB connection formulas around the classical turning point z = zb determine the ampli- 
tudes of region 777 as 



a 



U) 



iC, 



U) 



2Ti 



-i-n/4 



c y 7 



U) 



(j) 



2Vi 



+ C^Yi I e 



iir/4 



(3.22) 



where we introduced the quantity 



XnYj = Pi j \J\Qj\ dx 

J Z A 



(3.23) 
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At this point, we impose the second boundary condition of our problem. At the horizon, 
regularity of the perturbations in the ingoing Eddington-Finkelstein coordinates allows only 
ingoing modes. This requires setting C7 = in ( 3.22[ ). 

0) we find 



From the relations (3.17), (3.20), and (3.22) (with the boundary condition Cj 



that quasinormal modes must obey the condition 



liAT)-!)^^ + {AT) + l)e 



l , u 







& e~ iVj +ie ir >i ~0, for Pj > 1 , (3.24) 



where r)j and Tj are defined in (3.16) and (3.23), respectively, and in the last equality we used 
that ^> 1 in the large pj limit of our WKB approximation — see (3.23) 10 The quasinormal 



modes are thus described by solutions of (3.24). That is, the real part of the frequencies, 



of the quasinormal mode spectrum are quantized by the condition 



rjj(wj) =nir + 



7T 



Pj 



71 



W 



J, n 



Vdx= -(2n + 1 + Vj) 



n 



0,1,2, 



(3.25) 



with the non-negative integer n being the overtone that fixes the number of radial zeros of 
the mode, and the factor Vj that depends on the specific perturbation sector is defined in 



(3.14). At this point we just need to find the damping timescale of our system, as determined 
by the imaginary part of the quasinormal mode frequencies. Dissipation occurs because the 
bound states of region I can leak through the potential barrier of region II and be absorbed 
by the black hole horizon. This is thus a standard WKB tunneling process where the damping 
timescale Tj = (Imwj)" 1 is proportional to the inverse of the WKB transmission coefficient 
through the barrier, Tj ~ Tp i.e. 



exp ( 2pj 



~A 



V — Wj n dx 



n 



0,1,2, 



(3.26) 



where a prefactor can be explicitly computed but is irrelevant for our proposes since it is a 
finite, order 1, value. The massive scalar quasinormal mode spectrum in a large Is WKB 
approximation was first computed in [T3] (using an equivalent anti-Stoke line WKB formalism) 
and this WKB expression has been explicitly compared against a numerical computation in 
[T4"] . For large angular momentum 
large 



tj, the damping timescale (3.26) scales as Tj >n 
tj quasinormal modes are thus very long-lived. 



e l *. The 



We can solve (3.25) and (3.26) to get valuable physical information. As a first approxi- 



mation, for frequencies Wj L ~ 1, the bound states are localized near the boundary and we 



can approximate the potential well around z = by a simple parabola V 



L 2 



1 + 



with 



turning point at za — L\ w 2 -L 2 — 1. In this regime the analysis is essentially blind to the black 



hole presence and captures only the global AdS asymptotics. Inserting these relations into 
(3.25) we find that in this approximation the real part of the frequency, u n = Wj, n pj, naturally 



gives the normal mode frequencies of global AdS^, for the three sectors of gravitational modes 
and for the massive scalar field, 



Lw j>n ~ 



t , ,AdS 

t , .AdS 

t , .AdS 
LUJ S,n 
t , ,AdS 



(d-l) + 
(d-2) + 
{d-3) + 
A± + £ s 4 



7 T + 2n, 
W + 2n , 
7 s + 2n, 
2n, 



for tensor gravitational modes , 
for vector gravitational modes , 
for scalar gravitational modes , 
for massive scalar modes , 



(3.27) 



10 Large corresponds to a small imaginary part of the frequency, and in the end of our computation we 



confirm that indeed one has T'j ^> 1: sec (3.26). 
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Figure 3: Left Panel: Real part of the WKB quasinormal mode frequencies of the 
Schwarzschild-AdS black hole with respect to the AdS normal mode frequencies as a function of 
the WKB parameter Ej, in d = 4 (and overtone n = 0, and horizon radius r + = 0.1L). Moving 
from left/top to right/bottom, the curves describe the scalar gravitational (S), scalar field with 
BF mass fi 2 L 2 = —9/4, vector gravitational (V) and massless scalar field (/x = 0) cases (there are 
no regular tensor modes in d — 4) . Right Panel: Evolution of L Aus, n — L u>s, n~ L cu^^f as the 
dimension d increases, for the scalar gravitational (S) case. From the left/top to right/bottom 
of the figure we have the lines corresponding to the cases d — 6, d — 5, and d — 4. In these 
plots the red continuous line is the best fit curve of the data to the curve Auj tTl = ai?. The 
best fit value of /3 is — (d — 3)/2 to machine precision. 



for n — 0, 1, 2, • • • . In the last case, the choice of A + or A_ (defined in (3.9)) depends on the 
boundary condition for the scalar field. These expressions for the normal mode frequencies of 
AdS are known to be exact, i.e. valid for any value of angular momentum £j, in spite of the 
approximations used to obtain them. 

Beyond this global AdS (parabola) approximation, the real and imaginary parts of the 
quasinormal mode frequencies can be computed by numerically solving (3.25) and (3.26). We 
are particularly interested in the difference AcOj >n = uj^ n — w^jf. We find that in this large 
£j WKB limit, the real part of the Schwarzschild-AdS quasinormal mode frequencies has the 
universal power law behavior 



o>S,n~utt* + ^ j) £] 2 , (3-28) 

where is a constant that depends on the kind of perturbation (and radial overtone), 
dimension d and black hole radius in AdS units. The exponent — (d — 3)/2 was obtained by 
fitting Aujj tn to a power law for large £j. In particular, fitting to 1000 > £ > 600 for the values 
used in figure [3] yields —(d — 3)/2 to machine precision. Typical examples of Aojj )n data and 
the best fit curve are shown in the left panel of Fig. [3] for d = 4, overtone n = 0, and horizon 
radius r + = 0.1L, for the several perturbation sectors we study. On the other hand, in the 
right panel of this figure, we focus on the scalar gravitational sector of perturbations and study 
the evolution of the spectrum with the spacetime dimension d. A single illustrative example is 
shown, but we have explicitly checked the power law decay (3.28) for all perturbation sectors 
up to d — 11. When reading these plots one should bear in mind that the WKB approximation 



d-3 
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expressions are valid only for large £j, although numerical computation [2] shows the WKB 
expression for the scalar field to be a good approximation even for relatively small values of is- 



3.2 Long-lived quasinormal modes: Rotating black hole 



It is natural to ask what effects arise from introducing angular momenta, both in the background 
and in the azimuthal mode number. The most natural background in which to address this 
question would be the Kerr-AdS black hole or its higher-dimensional generalization, the single 
spinning Myers- Perry— AdS black hole. Unfortunately, an exact harmonic decomposition of 
general perturbations is not available for these backgrounds. So we instead choose to study 
the equal angular momenta Myers-Perry— AdS (MP-AdS) black hole in odd dimensions. As 
a consequence of its enhanced symmetry (it is cohomegeneity-1), exact analytic results are 
available both for the angular eigenfunctions and for the associated eigenvalues. This allows 
much of the WKB quasinormal mode analysis to be done analytically, with numerics playing a 
very minor role. The hope is that the key features of their spectrum, such as the deviations of 
the frequencies from AdS, extend to other rotating AdS black holes. 

The equal angular momenta MP-AdS black holes have the following line element [231 Wi 



where 



ds 2 = -f{r) 2 dt 2 + g{r) 2 dr 2 + h{r) 2 [d^J + A a dx a - Q{r)dt} 2 + r 2 g ab dx a dx b (3.29) 



9{r) 



r 

1? 



2M 

Z2N 



a 

L 2 



2Ma 2 



2-I- 1 



f(r) 



r 



g(r)h(r) 



n(r) 



r 

2Ma 



2N+2 



r 



2N h 2 



N 



h(r) 2 
rf-3 



2Ma 2 

r 2N+2 



(3.30) 



and g ab is the Fubini-Study metric on CP^ with Ricci tensor R ab = 2(N+l)g ab , and A = A a dx a 
is related to the Kahler form J by dA = 2J. 

Asymptotically, the solution approaches the Einstein static universe. The event horizon is 
located at r = r + (the largest real root of g~ 2 ), which is a Killing horizon generated by the null 
tangent vector d t + ftac^, where the angular velocity is 



n 



2Ma 



H 



„27V+2 



+ 2Ma 2 



< O ext 



where 



(~)Cxt 



1 + 



NL 2 



(N + l)r 2 



(3.31) 



The solution saturating the bound on Qh is an extreme black hole with a regular, but degen- 
erate, horizon. Note that the upper bound in Qh L is always greater than one and tends to 
unity in the limit of large r + /L. 



Consider now the Klein-Gordon equation for a massive scalar field (3.4) in this background, 



which was studied in detail in Assume the separation ansatz for the scalar field, 



r, x) = e- lujt e im ^ Y tm (x) h~* r~ N $ 



(3.32) 



L m 



where m must be an integer since ip has period 27r, and Y^ m (x) accounts for the dependence 
of the perturbation on the CP^ base space coordinates x. A remarkable (and very useful) 
property of the equal angular momenta MP black hole is that all the information about the 
eigenfunctions 3^ m (x) and its associated eigenvalues are known explicitly [271 [26]. This situation 
is to be contrasted with the Kerr-AdS or single spinning MP cases where the AdS spheroidal 
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harmonics and eigenvalues have no useful analytic expression. More concretely, l^ m (x) is an 
eigenfunction of the charged scalar Laplacian on CP^ [271 126] , 

(3.33) 



- V Y im (x) = A^ m F tot (x) , where V a = D a - imA a 
is the gauge-covariant derivative on CF N . Its eigenvalues are [27J [26] 



\ im = £(£ + 2N) - m 2 



2k + Iml 



fc = 0,l,2,... , 



Iml < £ 



(3.34) 



With the ansatz (3.32), the Klein-Gordon equation reduces to the angular equation (3.33) and 



to a radial equation for ^/ u im( r )- The latter can be written in the form of a time-independent 
Schrodinger equation, 



uiim 



(u - U+) (u - U-) $ 



ui£m 



o, 



Va 



PVh d 

r N+l dr 



f 2 h d 

r dr 



Vhr N 



(3.35) 



U± = mVL ± J V + 



P 



r 2 h? 



m?r 2 



+ h 2 



+ 2N) -m 2 + r 2 /i 2 + 4(1 - a) 



where a = 1 (for the scalar field; see below), and we have introduced the tortoise coordinate 



dr 



f 



< z < +oo 



(3.36) 



that describes the horizon at z — > +oo and the asymptotic boundary at z ~ L 2 /r — > 0. 

So far we have discussed a massive scalar field in the background of our rotating black hole. 
We would also like to consider perturbations of the background gravitational field. In order 
to disentangle perturbations of the Schwarzschild black hole, KI classified all perturbations 
according to how they transformed under isometries of the (d — 2)-sphere. Similarly, gravita- 



tional perturbations of the MP-AdS black hole (3.29) can be decomposed into scalar, vector 
and tensor parts, according to how they transform under isometries of the CP^ base space. 
The details of these perturbations are discussed in [26]. The relevant information for our own 
study is that, for tensor gravitational perturbations, solving the linearized Einstein equations 
reduces to solving a radial equation for a single gauge invariant quantity. This equation can be 
written exactly in the Schrodinger form (3.35) for /x = (and o ^ 0) [26j. Tensor perturbations 
exist only for odd d > 7, i.e. integer N >2. 



To summarize, the radial Schrodinger equation (3.35) describes perturbations of a massive 
scalar field (for a — 1 and d > 5) and, in addition, tensor gravitational perturbations (for /i = 
and d > 7)0 

We are interested in the large angular momentum WKB limit, £ — > oo, of the radial 



Schrodinger equation (3.35) . An important difference with respect to the non-rotating case 



discussed previously is that the behavior now depends on the azimuthal quantum number m. 
The most interesting WKB modes have \m\ ~ £ — > oo which we henceforward assume to be the 
case. Following closely the non-rotating computation of the previous section, we define 



^m(z) = e pS ^ z \ u = pw, £ = p-N, m = e^(e^ = ±l), Vj = ^/{N + l) 2 + fi 2 L 2 , 

(3.37) 



11 We refer the reader to [35] for a discussion of further constraints on a and on the eigenvalues that are 
irrelevant here. 
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where \ij = fi for the massive scalar field and \ij = for tensor gravitational perturbations. 
Equation (3.35) then reads 



P 



S"(z) + S'{zf - [V + (z) - w] [Vl(z) - w) 



1 

p 2 



Xjl 



with V±(z) = e^Q ± 



gh 2 



(3.38) 



and Xj{ z ) is a function associated to the higher order WKB contribution whose explicit expres- 
sion can be read from (3.35) but is irrelevant for our study. Notice that V± — > \{V defined in 



(3.11) when Qh — >■ 0. Note also that the harmonic decomposition of this subsection is done 

base space. When sending £l H — > and comparing results, we should 
N harmonics form a subset of the KI (S N ) harmonics. 



with respect to the 
have in mind that the 



The quasinormal mode frequencies in the limit of large p can be found by solving (3.38) 



in a WKB approximation. The procedure is now very similar to the Schwarzschild-AdS case 
analyzed in the previous section. The form of the potentials V± is illustrated in Figure |4j It is 
well-known that black holes with L > 1 are superradiantly unstable (i.e. its "quasinormal" 
modes grow in time in this regime). Here, we are interested in the case where this instability 
is absent, Qh L < 1; this is the case addressed in Figure |4| Under this condition, long-lived 
quasinormal modes do exist. They are bound states in the exterior allowed region / that 
can decay through the potential barrier (region II) downwards toward the black hole horizon 
located in the WKB region In the left panel plot, the upper black solid line describes the 
WKB potential V + (z) when the wave co-rotates with the black hole, i.e. when m — £ ^> 1. We 
compare it with the same potential in the limit Q — > 0, where we recover the Schwarzschild- 
AdS case of the previous subsection. This is the red dashed line. In the rotating geometry, 
the co-rotating modes see a potential barrier that is both higher and wider than the non- 
rotating one (for a given frequency w and £). Therefore, we anticipate that in the rotating 
background corotating quasinormal modes can be even more long lived than in Schwarzschild- 
AdS. The plot for V + (z) in the right panel shows that the opposite conclusion applies to counter- 
rotating modes that have — m — t ^> 1. To confirm and quantify these expectations we do 
the actual WKB analysis. Fortunately, we can borrow most of our calculations of the previous 



subsection, namely equations (3.13) and (3.15)-(3.24), as long as we take the replacement 



Qj(z) = to? — V(z) —> [vjj — V+(z)] [wj — V-.(z)], and Vj as defined in (3.37). 

In the end, we find the real part of the frequencies of the quasinormal mode spectrum are 
quantized by the condition 



7T 



7](Wj) = U 7T + — <^ p 



71 



Wj, n - V + ) (w j>n - V-) dx = -(2n + 1 + Vj) , n = 0, 1, 2, • • • 

(3.39) 

with Vj defined in (3.37) (fij = /i for the massive scalar field and fij = for the tensor 
gravitational perturbation). The damping timescale scales as 



T j,n ~ 



exp 2p 



~B 



~A 



w 



3,n. 



(V- 



w 



dx 



n 



0,1,2, 



(3.40) 



It follows that at large angular momentum the modes are again long-lived with a damping 
timescale (3.40) that scales as r ~ e e . Furthermore, evaluating (3.39) numerically and fitting 
Aojj tn to a power law again yields (3.28), though with constants which now depend on 
the angular momentum of the black hole background. We have checked this for d = 5,7, 9, 11 
and for both co-rotating and counter-rotating modes. This supports the conjecture that this 
behavior (with appropriately solution-dependent K d ) is universal for solutions asymptotic to 
global AdS. 
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m = t» 1 



-m = t» 1 




Figure 4: The black solid line describes the WKB potentials V±(z) for the rotating MP black 
hole (in the non-superradiant regime, Qh L < 1). For comparison, the red dashed lines describe 
the same potential when the rotation vanishes, i.e. the Schwarzschild-AdS limit. We also show 
the WKB allowed (/, and forbidden (II) regions for a given (rescaled) frequency w. The 
left panel describes corotating modes (m — £ ^> 1), while the right panel is for counter-rotating 
modes (— m — £ 1) 

4 Nonlinear stability of geons and boson stars 

We now investigate the stability of horizon-free stationary solutions such as geons and boson 
stars which are asymptotic to global AdS solutions. The modes most likely to drive such an 
instability are those with large angular momentum I and small radial excitation, since they 
have frequencies close to modes of AdS. Indeed, in the previous section we saw that for AdS^ 
Schwarzschild, such modes have frequencies ( |3.27| 3.28) 



d-3 , 



Lio = Loo AdS + 0(£-—)=£ + C + 0(£-—) (4.1) 

where C is a constant that depends on dimension and the type of perturbation. By examining 
a class of higher dimensional rotating black holes, we saw that rotation does not change the 
last term. This suggests that we will find similar behavior for perturbations about generic 
stationary solutions asymptotic to global AdS (see again footnote [8]) . In particular, horizon- 
free stationary solutions should exhibit good approximations to the resonances that drive the 
nonlinear instability in AdS. 

Whether or not there is an instability will depend on how much power is initially contained in 
these large £ modes. This is directly related to the degree of differentiability of the perturbations 
we choose to consider, and can be quantified as follows. Suppose we have a field on with 
Fourier coefficients 4>(k) ~ k~ p . To be in the Sobolev space H s we need 



d d ~ l k[k s (j)(k)f < oo. (4.2) 



This requires 



d — 1 , 
p>s-\ — . (4.3) 
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We certainly want to choose s so that our perturbations would induce no instability of Minkowski 
spacetime. In Einstein-Hilbert gravity, four dimensional Minkowski spacetime has been shown 
[28J to be stable for s = 3 and is further believed |29j to be stable even for s > 5/2. In d 



dimensions, it is believed [29] to be stable for s > (d+l)/2. From (4.3), this requires p > d. By 
the Sobolev embedding theorem, a field is in C q (i. e. q times differentiable) if p > q + d — 1. So 
the minimum requirement p > d ensures that the fields are in C 1 but not necessarily C 2 . We 
will find that geons and boson stars are nonlinearly stable unless they are in a high dimension 
and perturbations have low differentiability. 

We wish to add a small perturbation to a given solution and study its time evolution. 
We will address this problem by attempting to construct the general nonlinear solution using 
formal classical perturbation theory. By formal, we mean that we will focus on the existence 
of a good asymptotic series rather than attempting to demonstrate convergence. We will see 
that this level of argument predicts results that agree qualitatively with rigorous theorems for 
the nonlinear Schrodinger equation as reviewed in section [2j 

Recall that the perturbative expansion can be represented as a sum over connected tree 
graphs, where the vertices of the graph are associated with the monomials that appear in a 
Taylor series expansion of the action about the given background and the lines of the graph are 
associated with a given propagator. As usual, it is simplest to work in terms of normal modes. 
We assume that we are given a "seed" solution of the linearized equations of motion written as 
a sum over normal modes ip?, where I represents a complete set of labels for the normal modes: 
<f}W(x,i) = Yl/i 9r Yt We seek a solution of the form 



with 
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(f>{x,t) = Y^M t W x ^ ( 4 - 4 ) 



^) = E4 n) w- ( 4 - 5 ) 



n=l 



Here each (jy~' is an order n homogeneous polynomial in the first order coefficients (j>g and at 
each order n we require the equations to be satisfied only up to higher order terms. 

For n > 0, the terms of the polynomial (j)^ are associated with tree graphs with n + 1 

external lines as follows. For each such graph, label one external line with the pair (£, +). This 
line is to be thought of as the root of the tree (taken to be at the bottom) and is associated 
with the mode we are trying to compute. The other external lines are assigned arbitrary pairs 
(£j, ±j) for j = 1, . . . , n of modes and sign We may think of these as the leaves of the tree 
(and we take them to be at the top). Internal lines are also assigned arbitrary pairs of modes 
and signs. Each vertex v contributes a coupling constant g, an overlap factor C({£}), and an 
energy denominator 

v 1 -. , (4-6) 

zl^lines iEv ' 

where we say that % G v for a line % and vertex v if % is attached to v. The sum in the energy 
denominator has one term for each line connected to the vertex. For an external line or for 
the internal line below the designated vertex we set cDj = ijCUj where ±j, Ui are the signs and 



12 The n in this section should not be confused with the one in the previous section. Here n will represent the 
order of the perturbation expansion and not the radial excitation of the linearized mode. 

13 The inclusion of both signs stems from the fact that our modes are real and thus contain both positive and 
negative frequencies. 
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frequency associated with the mode assigned to that line. For an internal line above the given 
vertex v, we define u)j to be the (signed) sum ^ leaves j >ir over all external lines j that lie 

above % in the tree (j > i). In a theory without derivative couplings the overlap factor C({£}) is 
just an integral over space of the product of all modes associated with the lines connected to the 
given vertex, though in a gravitational theory two of these modes will appear with derivatives. 
These factors are multiplied together, and also multiplied by <jyu' . . . exp(— % YTj=\ 



The result is then summed over all non-trivia sets of pairs of signs and modes. For example, 



the second order term in any theory comes just from the cubic vertex and takes the form 

±i/i,± 2 / 2 



We now wish to estimate the result of the sum (4.5). As we will argue below that for d > 7 
the sum is dominated by the modes with no radial excitation, it is convenient to restrict to 
both d > 7 and no radial excitations from the beginning. At the end of our analysis we will 
also explain the result for 6 > d > 3. 

As mentioned above, we will focus on whether the sum associated with a given graph is 
finite; i.e. whether it is UV-convergent. Since we are interested in the high-momentum behavior 
we may use the flat-space approximation in which the sums over modes are replaced by integrals 
over Fourier space. This has the advantage that the overlap functions C({£}) become delta- 
functions that enforce conservation of spatial momentum, perhaps multiplied by certain powers 
of momentum associated with the derivatives in the given interaction vertex. Since all diagrams 
are tree graphs, these delta-functions can be used to eliminate the sum over internal momenta 
as well as one of the sums over external momenta. Since we study the sum with only modes 
on the sphere, each remaining momentum lives in R d_2 . 

It is useful to organize the momentum integrals in terms of an overall momentum scale k 
(defined, e.g., by k 2 = Y^=i l^«lV n ) an d a remaining set of integration variables {/%} which 
take values in g^- 2 )^- 1 )- 1 . The contribution of a given diagram of order n in the linear seeds 
(p^p can now be written 

graph„ ~ £ [ k^^dk [ dJfl^e-^A ]J *^4_, (4.8) 

sl^s J y S (d- a) (n-l)-l \ + ± *i / vertices ^Hnes i6v^ 

where % G v for a line i and vertex v means that i is attached to v. 

For Einstein-Hilbert gravity coupled to a scalar field, the largest interactions in the UV 
have two derivatives. So we may write F v (k, n) = k 2 F v (K,). We may also write (jfp = k~ p cf)i(K). 
Now, the gravitational action contains vertices of arbitrary order. However, as noted above, each 
vertex is associated with precisely two derivatives. On the other hand, we may use the cubic 
vertex to build an effective higher-order vertex of any order; figure [5] shows a 4-point example. 
Since the effective m-point vertex is built from m — 2 cubic vertices, it includes a factor of 
£2(m-2) f rom derivative couplings. While it also contains m — 3 extra energy denominators, 
each denominator is no larger than O(k). So the composite m-point vertex built from cubic 
vertices is at least 0(k m ~~ 3 ) times the size of the fundamental m-point vertex at large k. Thus 



14 A set is trivial if the corresponding sum vanishes no matter what frequencies are assigned to each mode. 
E.g. LUj; — ujj; + u!jr — uv at third order. Excluding these terms is associated with a resummation of naive 
perturbation theory which shifts the frequencies of each mode at each order but which should not affect the 



asymptotics (4.1 1 in the limit of small perturbations. 
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two 3-point 
interactions 



Figure 5: Two 3-point interactions combine to make an effective 4-point vertex. 



we see that the dominant graphs at large k involve only cubic vertices. Note that at order n in 



the linear seeds (jy- such graphs involve n — 1 cubic vertices We therefore neglect all other 
graphs for the rest of this section and write 

graph n ~ [ k^^dk [ dn (flhWe-vA ]J 9 ^ K \ . (4.9) 

signs J y S ( rf -2)Cn-l)-l \ j=1 J uev ^ices Klines iev "i 



It remains to estimate the energy denominators. From (4.1), we see that at large £ each 



Ui = O(k). A typical energy denominator thus also satisfies ^] lines iev Wj ~ 0(k). However, in 



some cases the 0(k) terms can cancel in the energy denominator. From (4.1) we see that this 
happens precisely when the sum Klines iev^* vanishes in the approximation Ui = £i, where £i 
is the integer associated with the total angular momentum of the mode i. The condition 
that this occurs may be used to express e.g. £2 in terms of the other angular momenta, 
reducing the integral over k to (d — 2)(n — 1) — 2 dimensions, and similarly for each additional 
constraint imposed. We note that exact resonances of this form were demonstrated in [H [3] for 
perturbations about empty AdS. 



An important observation is that, when this condition is imposed, the 0(1) term in (4.1) 
cancels whenever (when written directly in terms of external momenta!) the set of signs {±j} 
contains equal numbers of positive and negative signs. Note that making this choice of signs 
imposes no further constraints on the momenta. So, by making this choice of signs and imposing 
the above-mentioned single constraint we now find that Klines iev^« = 0(k~( d ~ 3 >' 2 ) or smaller. 
For this reason, momenta satisfying this constraint should dominate the result: While imposing 
the constraint reduces the integral by one dimension, this will "cost" no more than a factor of 
roughly k. On the other hand, we have increased the contribution from the energy denominator 
by k^ 1 ^ 2 . This is precisely the effect of the approximate resonances. 

Note, however, that such a choice of signs can be made only for energy denominators that 
involve an even number of external momenta. One might think that, when an odd number of 
external lines are involved, one could simply shift the constraint to require e.g. ^(^K* = C 
instead of = 0. However, the situation is not so simple. Indeed, refs. [U [3] studied 

the explicit form of perturbation theory about empty global AdS4 and found no resonances 
of this form at second order (three external lines). In simple cases (e.g. for three scalar-type 



15 At order n we have n + 1 external lines. Since each internal line connects to two vertices, connected tri-valent 
graphs satisfy E + 21 = 3V, where V is the number of vertices, E is the number of external lines, and / is the 
number of internal lines. But for any tree graph V = I + 1 by Euler's theorem. Thus / = n — 2 and V = n — 1. 
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gravitational modes) elementary considerations of SO(d — 2) representation theory force the 
associated overlap functions to vanish, removing any potential resonances. While we have not 
been able to rule out all potential 2nd order resonances in this way, we have explicitly studied 
all analogous cases involving a graviton (of any type) and two modes of a massive scalar field 
(of any mass) and found no resonances at this order. 

It is likely that new approximate resonances of the above form do in fact arise at higher even 
orders in perturbation theory. But these new resonances turn out to have little effect. We now 
make a short digression to explain why. The point is that, for a given diagram of any order, no 
more than half of the energy denominators can be approximately resonant in the above sense. 
This can be seen from the fact that our perturbation theory is dominated by diagrams built 
entirely from 3-point vertices as follows. 

First recall that, in computing the relevant u>i that contributes to an energy denominator, 
the definition of a), for a given internal line depends on whether this line lies above or below 
the vertex in question. However, for a given labeling of the external lines, the two definitions 
coincide precisely when vertex above the line has an exact resonance. When this coincidence 
occurs, we refer to the line in question as a resonant line. We may similarly speak of approxi- 
mately resonant internal lines using the above sense of approximate resonances. Now, imagine 
dividing a given diagram into subdiagrams by cutting all of the approximately resonant internal 
lines. The subdiagram containing the original root merits special treatment: if the root (an 
external line) is non-resonant we will discard this subdiagram, though we keep it if the root is 
resonant. 

With this understanding the (approximately) resonant energy denominators for the original 
diagram are in one-to-one correspondence with the subdiagrams obtained above, so that count- 
ing one also counts the other. Due to the above coincidence regarding &i, each subdiagram 
comes with a definite labeling of its external lines by (perhaps approximate) modes. Since by 
construction each such subdiagram contains a resonance, it must contain more than one 3-point 
vertex. That is, the number of sub-diagrams is no more than half of the number of vertices. 
This verifies the claim that the number of approximate resonances is also bounded by half the 
number of vertices, i.e., there are at most n of the above approximate resonances at either order 
2n + 1 (with 2n vertices) or 2n + 2 [2n + 1) vertices. 

Returning to the main calculation, the above argument motivates what we call the generic 



frequency assumption. Namely, we will assume that the integral (4.8) is dominatec 16 by con- 
figurations in which, for each of the (n — 1) vertices, the constraint ^j(±)i^ = is satisfied. 

We furthermore assume that the integral is dominated by configurations (of momenta and 
signs) for which the contribution due to each energy denominator with an even number of 
terms is 0(k( d ~ 3 " 2 ) and for which the contribution due to each energy denominator with an 
odd number of terms is 0(1) at large k. We then write Klines iev ^ = ^ v (k), k~^ d ~ 3 ^ 2 S v (n) 
as appropriate. As noted above, there are at most n approximately resonant vertices at either 
order 2n + 1 or 2n + 2. We therefore focus on odd orders (2n + 1) which will clearly show the 
strongest UV effects. 

Letting k denote configurations that satisfy the constraints associated with the above- 

16 Or, at least 'marginally dominated' in the sense that the contributions from such configurations are of 
the same order at large k as the full answer. That is, we may neglect other sets of configurations which give 
contributions of the same order, such as those obtained by not imposing the constraint XL(±i)-^j = f° r an 
energy denominator with an odd number of external momenta. When all phases are aligned, this would allow 
the phase space volume to contribute an additional factor of k and would reduce the contribution of the energy 
denominator by only a compensating factor of k. 
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mentioned generic frequency assumption, we may now estimate (4.9) at order 2n + 1 to be 



2n , 

3d-*- v 



gra P h 2n+1 ~ [ k n{d-S)/2 k 2n d - { 2n + l) {p+ l) dk f / \\^{R) e ^A TT 

J ,/ S (d-3)(2„-l)-l W / . <er ±£ 



Su(/c) 

DGvcrtices v ' 



(4.10) 

where T^ d 3 ^ 2n 1 is a (d — 3)2n — 1 dimensional manifold defined by imposing the required 
constraints on g^ -2 )^ 71-1 ) -1 . A similar result holds at order 2n. 

We are now in a position to justify our restriction to modes with no radial excitation for 
d > 7. A study of black hole quasinormal frequencies for modes with radial excitations indicates 
that they admit an expansion Aw = I + Cq + C\/l + . . . at large i with coefficients Cj that 
depend on the radial quantum number. Let us consider the effect of radial excitations at third 
order in perturbation theory, where we have one energy denominator that refers to 4 modes. 
As usual, we take two to be of each sign to enhance cancellations. However, the constant term 
Co cancels only when the radial quantum numbers agree pairwise. That is, the two 'positive' 
modes have the same radial quantum numbers as the two 'negative' modes. This imposes two 
constraints on the 3 sums over radial quantum number. So the extra phase space associated 
with the sum over radial modes with radial quantum number <C I can contribute no more than 
a factor of £. Furthermore, we generally find an energy denominator of order l/£ as opposed 
to the l/£~i~ obtained in the case with no radial excitations. For d > 7, we see that the sum 
is indeed (at least marginally) dominated by the restricted sum over modes with no radial 
excitations. The same effect occurs at higher orders. 

The last ingredient to consider are the phases e~ luJjt and the phases in <^(/c). The result is 
largest when all of the oscillatory terms are in phase; e.g. when the <j)je~ luJ:it are all real and 
positive. This may take an exponentially long time to achieve, but since we are interested in 
late time stability, we will estimate the result assuming this condition is satisfied. 

In this case the k integral becomes independent of k and we find 

graph 2n+1 ~ j k n{ - d -V/ 2 k 2nd - {2n+1){ - p+1) dk. (4.11) 
The integral over k converges in the UV if 

(d-3)/2 + 2(d- 1) -2p<p/n. (4.12) 

So for 

5d — 7 , 

P>^~ (4.13) 

it converges for all n while it diverges at large n for smaller p. The terms at order 2n are slightly 



more convergent at large k, so the UV properties are determined by the odd terms (4.11). Note 
that when the integral converges in the UV, terms of higher order in n scale with higher powers 
of the seed solution's amplitude and we have a good asymptotic series. In such cases it then 
appears that the solutions are nonlinearly stable over arbitrarily long time periods for initial 
data that lies in a Sobolev space H s for s > (3d — 5)/4 (using ( 4.3[ )). The similarity of this 



conclusion to theorems regarding the nonlinear Schrodinger equation (reviewed in section [2 
encourages the belief that it can also established more rigorously. 

Recall that to avoid an instability in Minkowski spacetime, we need at least p > d. Our 
analysis above indicates a divergence with p > d when (hd — 7)/4 > d or d > 7, i.e. starting 
at AdSs. But as discussed at the beginning of this section, even that requires perturbations 
which are only C l and not C 2 . To see an instability with perturbations that are C 2 , one needs 
to go to at least to AdSi 2 . One may in fact need to go even higher. As noted in [3], the 
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detailed structure of the gravitational interactions can cause certain overlap coefficients C({£}) 
to vanish, removing potential (in this case approximate) resonances and making the sum more 
convergent than our analysis would indicate. It is also possible that further resummation of the 
series (analogous in the case of exact resonances to introducing additional frequency shifts at 
each order in perturbation theory) might further improve convergence. One should thus take 
our results as giving a lower bound on the dimension required for instability at a given degree 
of differentiability. 

The analysis above was valid only for d > 7. For smaller dimensions one must include the 
sum over radial quantum numbers. As noted above, this contributes an extra factor of k~ 
at third order. For each additional two orders, it contributes an extra factor of k from the 



d-S 



additional phase space, but reduces the integrand by a factor of k 2 due to the larger energy 

n(T-d) ■ ■ 

denominator, i.e., at order 2n + 1 it contributes k 2 . The analog of (4.13) is then p > d, 
which is the same condition under which Minkowski space is expected to be unstable. 

Even in cases of high dimension and low differentiability, one needs to wait for the phases to 
line up to see an instability. A more detailed estimate of (4.10) which includes random phases 
shows that the integral over k always converges in all dimensions assuming only the minimum 
differentiability required for stability of Minkowski spacetime. 

For AdS black holes, since the large i modes are long lived, we can ask if the above phases 
are indeed likely to line up before the relevant modes decay. Since the probability of N phases 
lining up to accuracy e is of order e N , we will generically need to wait a time of order e~ N 
before the phases line up. Since there are /c 2n ( d_3 ) phases in (4.10) at momentum scale k, these 
phases will line up only at ti n phase ~ t ~ k2n d 3 . But on the other hand we have seen that in the 
presence of a black hole the decay time for modes of momentum k is only e k <C t- m phase- So for 
d > 4 the phases are not in fact likely to line up before the relevant modes decay. 



5 Non-coalescing binaries 

In this section we describe a new class of asymptotically AdS solutions consisting of non- 
coalescing black hole binaries. It is easy to see that such solutions should exist. For simplicity, 
we consider four dimensional spacetimes. Consider a non-extreme Kerr AdS black hole with 
angular velocity Q > l/L. The Killing field £ = d/dt + Qd/dtp which is null at the horizon 
is timelike just outside the horizon and spacelike near infinity. So there must be a radius at 
which A = evaluated on the equatorial plane, has an extremum. The 9 — >■ n — 6 symmetry 
ensures that at this radius, V^A = 0. Now 

£ M V M r = -^A = (5.1) 

at this radius, so the Killing field is actually tangent to a geodesic. This geodesic represents a 
circular orbit around the central black hole. Imagine placing a small black hole on this circular 
orbit. This will perturb the background geometry, but since the orbit of the small black hole 
is invariant under the Killing field £, the perturbation generated will also be invariant under 
this Killing field. Since £ is the horizon generator of the original Kerr AdS black hole, this 
means that no energy falls into the black hole at the center. The radiation produced by the 
small black hole forms a standing wave which just supports the orbit. We have constructed 
this radiation explicitly and it will be presented elsewhere [30J. It should be possible to add 
higher order corrections to this perturbative argument and find exact solutions of this type. 

One can ask whether the above circular orbit is stable to radial perturbations. It turns out 
that if the black hole is too close to extremality, the orbit is unstable, but there is a range 
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of parameters for Kerr AdS so that the orbit is stable. Unfortunately, even when the orbit 
is stable, the entire solution is still unstable due to superradiance. This is because we need 
Q > 1/L to have a circular orbit which is invariant under the horizon Killing field. Generic 
perturbations will probably cause the central black hole to spin down and the orbiting black 
hole will spiral in as usual. 

One might also ask if one can place a small black hole in permanent orbit around a geon, 
where the radiation produced would again support the orbit against spiraling in. The answer 
appears to be no. The Killing field of the geon, £ = d/dt+ (uj/m)d/d<p is timelike at the origin 
and spacelike at infinity. But since there is no analog of the horizon where the Killing field is 
null, A = need not have any extreme 17 

It is amusing to note that for the non-extreme Kerr metric (in asymptotically flat spacetime) 
the horizon generator £ = d/dt + Qd/dip is always timelike near the horizon and spacelike near 
infinity for any Q ^ 0. So there is always a circular geodesic which is invariant under the Killing 
field, and putting a black hole on this orbit could generate radiation which is invariant under 
£. But the corresponding solution is not asymptotically flat as the radiation both to and from 
infinity is periodic in time. Instead of generating this solution, starting the binary on a similar 
orbit with asymptotically flat boundary conditions must lead to coalescence. Clearly, the AdS 
boundary conditions play a crucial role in the above construction. 



6 Discussion 

We have argued that many solutions asymptotic to global AdS are nonlinearly stable, despite 
the fact that global AdS itself is unstable. The key difference is that normal modes in global 
AdS have frequencies which are all integer multiples of the AdS frequency, so there are infinite 
towers of resonances in the perturbation expansion. In contrast, other solutions with less 
symmetry do not have exact resonances. While approximate resonances (associated with high 
angular momentum modes) still exist, we have seen that they are not strong enough to cause 
an instability. Our arguments were based on formal perturbation theory, though the qualitative 
similarity between our results and established theorems for the nonlinear Schrodinger equation 
encourages the belief that they can be placed on a more rigorous footing. There is also numerical 
support for the nonlinear stability of boson stars [3T] • The suggestion [3] that there may be a 
theorem showing generic solutions asymptotic to global AdS are singular is clearly incorrect^ 
While our methods allow us to study only the stability of stationary backgrounds, we expect 
perturbations of oscillon solutions to behave similarly. 

Particular solutions may of course have excitations with additional exact or approximate 
resonances at finite frequency. But as discussed in section |2j nonlinear instabilities typically 
require an infinite tower of (perhaps approximate) resonances. Since excitations of backgrounds 
asymptotic to global AdS have a discrete spectrum, any such tower involves modes of arbitrarily 
high frequency so that our analysis should suffice^ 

17 A continuity argument shows that there are no extrema for geons with small energy. 

18 See [35] for some of the difficulties in trying to prove such a theorem. 

19 Solutions asymptotic to Poincare AdS can also feature what is effectively an infinite tower of resonances 
at low frequency. In particular, at leading order in long wavelength along the AdS boundary, low frequency 
perturbations of planar AdS black holes are described |33j by the Euler equations for a conformal fluid on the 
boundary spacetime. These equations have exact zero modes associated with divergence-free fluid flows. It is 
natural to expect such zero modes to lead to nonlinear instabilities of planar AdS black holes. While black 
holes asymptotic to global AdS should be nonlinearly stable for small enough perturbations, fixing the size of 
the perturbation while taking the limit of a large such black hole should result in behavior much like that of 
the nonlinearly unstable planar case. We believe this explains, from the gravitational side, the results of [34] . 
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This argues that various "excited states" such as geons and boson stars can be stable. How- 
ever, one expects to also construct stable ground states by modifying the boundary conditions 
at infinity so the lowest energy solution differs from AdS. This includes putting a metric on the 
sphere at infinity which is not conformally flat. In theories with scalar fields whose mass is close 
to the Breitenlohner-Freedman bound, one has a large class of allowed boundary conditions for 
the scalar field. In many cases, the ground state includes a nonzero scalar field (35] . Finally, the 
AdS soliton [36] is believed to be the lowest energy solution for its boundary conditions. All of 
these solutions are likely to be nonlinearly stable by the arguments given earlier. In fact, they 
fall into the category of solutions mentioned in footnote [8] which should be even more stable 
than the geons, boson stars, and black holes addressed in the main text. 

The picture that emerges is that AdS is very special. Its high degree of symmetry leads to 
a large number of resonances and causes its nonlinear instability. Most other solutions do not 
have this feature. 

We now consider the implications of these results for gauge/gravity duality. It was argued 
in [3] that the instability of AdS has a simple interpretation in the dual field theory: Since 
black holes have the maximum entropy in the classical limit, the instability toward forming 
small black holes can be viewed as thermalization (in a microcanonical ensemble since the 
energy is fixed). The existence of geons was puzzling from the dual field theory viewpoint, 
since they represent finite energy excitations of the strongly coupled large N field theory which 
never thermalize. We now see that AdS is very special. Not only are geons stable, but so are 
most nontrivial ground states. This would seem to contradict the idea that energy added to 
a system generically thermalizes. Note that since we have finite backreaction in the bulk, the 
energy we add is always 0(N 2 ) and is not really small. So this lack of thermalization is not 
just a property of a few low lying energy eigenstates. However, one must remember that these 
results only hold at large N. In many ways, large N field theories act classically, and classical 
field theories almost never thermalize. If there is a potential landscape, a small amount of 
energy added near one local minimum will stay near that minimum and never explore other 
configurations. From this viewpoint, pure AdS is like having no potential. 

At finite N, one expects quantum tunneling to cause the states of the boundary theory (on 
S n x K) to explore all accessible regions, so many operators will behave thermally. Similarly, 
in the bulk we expect added energy to generically thermalize, (although the maximum entropy 
configuration can now be a gas of radiation and not a black hole). 
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